Asymptotics of the dispersion interaction: analytic benchmarks for van der Waals 

energy functionals 
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We show that the usual sum of R~ 6 contributions from elements separated by distance R can give 
qualitatively wrong results for the electromagnetically non-retarded van der Waals interaction be- 
tween non-overlapping bodies. This occurs for anisotropic nanostructures that have a zero electronic 
energy gap, such as metallic nanotubes or nanowires, and nano-layered systems including metals 
and graphene planes. In all these cases our analytic microscopic calculations give an interaction 
falling off with a power of separation different from the conventional value. We discuss implications 
for van der Waals energy functionals. The new nanotube interaction might be directly observable 
at sub-micron separations. 
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Dispersion interactions (part of the van der Waals, 
vdW energy) are especially significant in soft matter. 
The vdW physics that we expose here could be rele- 
vant in predicting the energetics of bundles of metal- 
lic nanowires or nanotubes, layered metallic systems, 
7T— conjugated systems including graphite, intercalated 
graphite, graphitic hydrogen storage systems and pi- 
stacked biomolecules, and other weakly bound ("soft") 
layered and striated nanosystems. Standard local (LDA) 
and gradient (GGA) density functionals for the electronic 
energy do not obtain any distant dispersion interaction, 
but density functionals have been derived recently that 
obtain, in a natural fashion, both distant dispersion in- 
teractions and their saturation at small distances. These 
and other numerically practicable vdW energy schemes 
available to date 0, M IE S IE IE 13 for the above 
systems (in the electromagnetically non-retarded regime) 
have a "universal" feature: the distant vdW interac- 
tion energy between sufficiently separated subsystems is 
given qualitatively by a sum of contributions of form R~ 6 
between microscopic elements separated by distance R. 
This leads to "standard" power laws E oc —D~ p for the 
interaction energy between various macroscopic bodies 
separated by distance D (column 3 of Table 1). Although 
these "universal" asymptotic results are indeed valid for 
most macroscopic systems, we show below that they fail 
for the anisotropic nanostructures mentioned above. Col- 
umn 2 of Table 1 summarizes the asymptotic (D — > oo) 
benchmarks that we propose below for the vdW energy 
of two parallel nanostructures of infinite extent. 

To analyze these situations, we use the correlation en- 
ergy E!} PA (D) from the Random Phase Approximation 
(RPA) [II EE III EE EE El, a basic microscopic theory 
that does not rely on assumptions of locality, additivity 
nor i? -6 contributions. Going beyond the RPA does not 
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TABLE I: Asymptotic vdW energy of parallel structures. K 
and Do are constants. * denotes new derivations given here. 



change the asymptotic power laws predicted here, un- 
less the exchange-correlation kernel f xc [lH 115 , Il9| has a 
slower spatial decay than the bare coulomb interaction, 
an unprecedented and unlikely scenario. 

Where the separated subsystems exhibit lightly 
damped long-wavelength plasmons, we note[2(| that the 
principal contribution to E^ p (D) comes from the sum 
of coupled- plasmon zero-point energies: otherwise we use 
the full RPA. Some essential common features of these 
systems will be abstracted from these specific calcula- 
tions. We obtain analytic results for the asymptotic 
(D — > oo) regime in all cases, but in section E we will 
also discuss systems near their equilibrium spacing. 

A: Distant attraction between metallic linear struc- 
tures. Consider two parallel, infinitely long conducting 
wires or tubes separated by a distance D substantially 
exceeding their radius b, and with b < A where A is a 
bulk screening length. Both standard ^ i?~ 6 analysis^ 
based on the vdW interaction between electrons local- 
ized in atoms or bonds, and recent functionals |2l| . give 
a vdW energy per unit length of the form E cx — D~ 5 . In- 
stead we consider the zero-point energy of the delocalized 
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coupled one-dimensional plasmon modes with wavenum- 
ber q parallel to the long axis 0,H2 .The radially-smeared 
intra-wire coulomb interaction is Wn(q) — w 22 (q) = 
— 2e 2 \n(qb) where 1 and 2 refer to the two wires, and 
we have assumed qb << 1, as appropriate when D >> b. 
In the same limit the bare density-density response for 
electronic motion parallel to the wire is Xon — X022 = 
N^q 2 1 \mu) 2 ) where wis the frequency, No is the number 
of electrons per unit length and m is the electron mass. 
RPA screening yields the interacting response of a single 
wire as xn = Xon/(l — wnXon)- The inter-wire coulomb 
interaction in the present limit has a Bessel form, w\ 2 = 
2e 2 Ko(qD) . The RPA equation for coupled ID plasmons 
on two identical wires is Xii w i 2 = 1; giving two roots 
for each q: w±(D) = c 1D \q\ (|ln(g6)| ± K^qD)) 1 ' 2 . Here 
c\d = (2iVoe 2 /m) 1 / 2 is a characteristic velocity. The 
vdW energy is the separation-dependent part of the sum 
of zero-point plasmon energies per unit length: 
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For D » b we expanded to 2nd order in 
Ko(qD)/ \ ln(qb)\ which is small near the peak of the in- 
tegrand. This gave |j| approximately 

E vdw /L « -(16^)~ 1 nci D I?- 2 (ln(2.39D/6))- 3/2 . (1) 

This approach is reasonable when the electron mean free 
path do al ong the wire satisfies do > D. In fact bismuth 
nanowires[23| and conducting nanotubes 24] can both 
have do > 1 micron. Eq. differs from the widely 
accepted result E oc —D~ 5 by nearly three powers of 
D, and is necessarily dominant, at sufficiently large D, 
over any such higher-power contributions (arising from 
the remaining bound sp 2 electrons (in nanotubes) and 
azimuthal 7r plasmons). Our plasmon model does give 
D~ 5 if a pinning force is added to mimic an insulator 0. 

B: Distant attraction between thin conducting layers 
Consider infinite parallel metallic plates separated by dis- 
tance D and of thickness b, with b « D, b < A where 
A is the bulk screening length. As is already well-known 
0,0], the zero point energy of long- wavelength coupled 
2D plasmons leads to an attraction of form E oc —D~ 5 / 2 . 
The ^2 R~ 6 approach, correct for thin insulators, gives 
E oc — £>~ 4 , different by 1.5 powers of D. 

C: Distant attraction between planar n-conjugated sys- 
tems. What does the physics of long- wavelength ex- 
citations imply for the energetics of layered planar n- 
conjugated systems, such as the controversial 
Wk Efj j and technologically important 0, E3 graphene- 
based systems? Firstly, an isolated graphene layer at 
T=0K is not a metal but a zero-gap insulator [3]. Thus 
one cannot argue for a metallic — Z)~ 5 / 2 energetics (as un- 
der (B) above) at large layer separation D and T = OK, 
even though band overlap makes graphite weakly metal- 
lic at the equilibrium layer spacing. We briefly derive be- 
low, however, our new result that the attractive energy 



between two well-separated graphene planes at T = OK is 
of form —C%D~ 3 , closer to metallic D~ 5 / 2 behavior than 
to insulating D~ 4 behavior. All the new physics here 
comes from electrons close to the Fermi level: we can ig- 
nore the response of the tightly-bound covalent sp 2 elec- 
trons, whose finite energy gap ensures that they produce 
a conventional vdW attraction of 2D insulator type (en- 
ergy oc — D~ A ), negligible at large separations compared 
with the D~ 3 vdW attraction that we shall find between 
the 7r 2 electrons of interest here. The bonding and anti- 
bonding 7r bands have a gapless bandstructure0|- The 
energy near the K points where the bands touch is given 
by e^ 1 ' 2 \p) — ^fhvo \p\ where p is the 2D crystal momen- 
tum measured from a K point, and vo is a characteristic 
velocity (about 5.7 x 10 5 m/s for graphene). From 
turbation theory within a Wannier description 
the zero-temperature density-density response Xks of in- 
dependent ir z electrons moving in the groundstate Kohn- 
Sham potential of a gapless 7r-layer is then of the form 
XKs(q, 0, 0, z, z', uj = iu) = S(q, z]S(qjZ')*Xo(q, iu), with 
/ Sdz — ► 1 as q — > 0. We found bU| the effective 2D 
response xo at small q and imaginary frequency ui = iu 
to be 



Xo(q, iu) w -2hv q (l + u 2 / {voq) 2 



-1/2 



(2) 



consistent with previous real-w results 32] . Here we treat 
the response in each sheet as strictly two-dimensional, 
and ignore certain local-field effects, so that the only 
consequence of the periodic potential is to replace the 
2D free-electron bare response —noq 2 /(mu 2 ) by the zero- 
gap Bloch response J5J). This is justified esewhere 
|3fl IsH ] . We consider electron density perturbations of 
form ni exp(iq.r + ut) in layer #1, where <fand r are 
two-dimensional. Such charge disturbances interact via 
a Fourier transformed bare coulomb potential, 

wn\(q) = 2ir\e 2 q~ 1 , w 12 \{q) = 2Tr\e 2 q~ 1 exp(-qD), 

(3) 

for interactions within a layer and between two layers 
distant D, respectively. Then the RPA equation for the 
interacting density fluctuation in layer #1 as driven by an 
external potential vf xt exp(iq.r+ut) is of time-dependent 
mean-field form, n± = xo(q, iu) (vf xt + wnni). This 
applies in the absence of layer #2 or equivalently for 
D — > 00. Solving for niwe find a single-layer density- 
density response 



Xiia,d^oo = ni/vl xt = xo/(l - wnxxo)- 



(4) 



With two layers present, the density response obeys cou- 
pled RPA equations n x = Xn\,D->oo(vf xt + w 12 \n 2 ), 

w 2 i\ni). The solution is 
(ni,n2) T and similarly for 
v ext , while the components of the 2x2 matrix x are 
Xm,£> = XiiA,r>->oo/(l ~ wi2aXiia,d^oo) and xi2\,D = 
w 12AXiia,_d^ooXiia,d- For the case of two identical layers 



n 2 = X22\,D^oo(v 2 xt 
n = xu ext where n 
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considered here, the other elements are X22X.D — Xnx,D 
and \2i\,D = Xi2\,D- 

In the present system the response © of a single layer, 
continued to the real frequency axis, yields no lightly 
damped plasmons (poles) for small q, so that a sum of 
plasmon zero-point energies cannot be used to evaluate 
the vdW interaction. Instead we consider the electro- 
magnetically non-retarded groundstate electronic corre- 
lation energy, which for a general inhomogeneous elec- 
tronic system is given exactly by the adiabatic connection 
fluctuation-dissipation theorem (see e.g. |1(J): 



E r = — 



2n 



dX / dfdf'- 



Axx(r,f',iu)du. 

(5) 

Here Axx = Xx ~~ Xo, where xx is the electron density- 
density response function at reduced coulomb interaction 
Ae 2 / \f — r'\. Applying © to the present layer geome- 
try and Fourier transforming parallel to the layers we 
find that the separation-dependent part E vdw /Aoi the 
energy per unit area is: 



E C (D) - E c {oo) 



A 



h - r d uf dx 

* Jo Jo 



2nqdq 



x {wnx (x 



11XD 



A Jo (2tt) 2 

XllA.D^oo) + Wi2\X12\d) • (6) 



Within the RPA approximation, Eqs. © and @, plus 
(0J and the equations following it, show that each term 
of form w%in 10 depends on u solely through the di- 
mensionless combination x = u/(voq). The remaining 
dependence of wx on q is solely via y = qD. Thus © has 
a scaling form 

fX i\ />oo />oo 

E vdw /A = h / a* / i d —,qD) 

Jo A J J v q 



too 
D 3 



1 dX 

T 



y l dy I dxG(X, x, y) (7) 
lo 



where G(A, x, y) is independent of D. We numeri- 
cally evaluated the dimensionless 3D integral in © for 

graphene parameters sua, 

giving the interaction en- 
ergy per unit area in Gaussian esu units: 



E vdW /A 



-7.745 7 xlO _2 to -D" 3 



-2.003 6 xl(T 2 e 2 Z)- 3 
(8) 

This D 3 form shows that the gapless 7r-conjugated 
planes behave in this respect more like metals (E oc 
— Z?~ 5 / 2 ) than insulators (E cx — D~ 4 ) , despite the lack 
of undamped 2D plasmon modes on a single n sheet. 
D. Parallel metallic and it-conjugated planes 
Another interesting case is the interaction between a 
7r-conjugated layer and a metallic 2D layer with fermi 
energy ep (e.g. an undoped and a doped graphene 
sheet). For D » D Q = h 2 v^/(2Tre 2 e F ) (=0{lnm) for 
ep = 0(0.02 eV) as in a bulk graphite layer) the methods 
described above give an energy per unit area (c.f. ©) 

E vdw j A _ -Ce 2 D- 3 ln(D/D ) (C constant). (9) 



As in the case of two non- metallic gapless ir layers, the 
result disagrees with standard theories. 

E. Interaction energy near overlap We now dis- 
cuss possible difficulties with the non-asymptotic, near- 
equilibrium energetics of the present systems, especially 
graphenes. The commonest ab initio approach, the LDA, 
misses distant dispersion interactions entirely |33l| . and 
yet gives a good lattice spacing pjif and good breath- 
ingphonon frequencies [35| in g ra p hite, (unlike GGAs 
|j,|3jj(]). Recent experiments |25|,|2jj , however, lead one 
to suspect jl^l that the LDA pays for its neglect of dis- 
persion physics by severely underestimating the equilib- 
rium binding energy of gra phite. LDA also underbinds 
related fullerene systems [27]. This phenomenon has been 
investigated in layered jellium analogs, via fully nonlocal 
many-body correlation theory (|13|i Fig. 4 of [13 )■ It 
was found that either layer-layer forces or binding en- 
ergy have serious errors near equilibrium, when distant 
dispersion forces are underestimated. The underestima- 
tion is related in turn to the lack of distant correlated 
fluctuations, especially those oriented parallel to the lay- 
ers. Thus these low-q fluctuations can have effects even 
near the equilibrium spacing. Addition of explicit R~ 6 
vdW terms has been a common remedy for stretched 
graphitic systems 0, 1 an d recently several seamless 
vdW schemes have been proposed 0, Q, 0] , based on 
approximations for response functions. Refs and Q 
are the mostpractical, and are qualitatively successful in 
graphitics. |j| overestimates the binding energy of small 
systems but correspondingly obtains a large binding en- 
ergy of two graphene layers (more than twice that from 
||| or from LDA, and consistent with experiment). Ide- 
ally a single theory should give reliable results for small 
and extended systems. Could it be that the key is a 
correct treatment of the fluctuations parallel to a long 
axis in the extended cases, the same fluctuations respon- 
sible for the unusual asymptotics exposed here that is ab- 
sent in 0,0? These fluctuations are of course dominant 
only at large separations but they might not be negligible 
near the equilibrium spacing, where all wavelengths can 
contribute. We speculate further that the same physics 
might apply in other large finite 7r-conjugated systems 
(e.g. planar melanin layers, carotenes, fullerenes |27j ) 
where, as the system size increases, the electronic gap 
diminishes while longer-wavelength excitations become 
possible. 

Summary and Discussion. Our new results (see 
0!©,© and Table 1) show that that usual sum of 
CqR~ 6 terms incorrectly predicts the dependence of the 
dispersion energy on separation D for a range of sys- 
tems. Simple energy functionals presently available all 
have standard C*6i? _6 asymptotics. A finite sum of 
multipole, or triplet and higher terms will also not re- 
produce what we have discussed. The standard asymp- 
totics fails when the component systems (i) are metal- 
lic (or have a zero electronic Bloch bandgap), and (ii) 
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are spatially extented in at least one dimension, so that 
long-wavelength (low-g) charge fluctuations can occur, 
and (iii) are of nanoscopic dimensions in another spatial 
direction, so that the electron-electron screening is re- 
duced compared with 3D bulk metallic systems, leaving 
a divergent screened polarizability at low frequency and 
wavenumber. (Thick metal slabs, for example, violate 
(iii): they have complete screening and exhibit a conven- 
tional power law, E oc —D~ 2 . See e.g. (12). Where 
free \ow-q plasmons exist, conditions (i) - (iii) imply that 
they will be gapless. The same conditions ensure that 
the usual spatially local approximation for the dielectric 
function [2j is invalid. Our results provide unequivocal 
asymptotic benchmarks that are not satisfied by exist- 
ing simplified van der Waals energy formulae, because 
they do not treat in enough detail the fluctuations along 
the extended space dimension. In Section E we have 
further motivated the possibility that the same fluctua- 
tion physics may be relevant in the systems considered 
here, even near their equilibrium spacing. Investigation 
of this question requires a seamless energy formalism that 
is fully nonlocal - e.g. RPA-like theories [Til [TH fl7l l20| . 
Such calculations are only now becoming possible for 3D 
systems [IE , with no converged results available to 
date for the present zero-gap cases. Simplified vdW en- 
ergy functionals are therefore certainly needed for routine 
modelling, and the above considerations suggest that ex- 
isting functionals may need further refinement to take ex- 
plicit account of large-scale geometry and/or nonlocal en- 
tities such as electronic bandgap|20|. We note finally that 
our work predicts novel differences in the forces between 
conducting and nonconducting nanotubes or wires, that 
might be directly measurable for low-index nanotubes at 
sub- micron separations 0, and that could even affect self- 
assembly processses. These considerations might also af- 
fect the analysis of some seminal experiments [jjl Efij con- 
cerning graphitic cohesion, because these relied at some 
point on theory involving a sum of R~ 6 contributions. 
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